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ABSTRACT 

[unclassified] 


i!:  a 


An  analytical  study  is  described  of  the  phenomena  leading  to  rod  loss 
in  the  perforation  of  finite  targets  by  hypervelocity  impact  of  end-oriented 
rods.  The  investigation  is  based  upon  two-dimensional  hydrodynamic  PIC 
(particle -in-cell)  code  solutions  of  sixteen  impact  cases.  These  impact 
cases  include  solid  and  hollow  rods,  normal  and  oblique  incidence,  and 
aluminum  and  iron  rods  impacting  aluminum  plates  of  varying  thicknesses 
at  velocities  of  4.6  and  6.1  km/sec  (15,000  and  20,000  ft/sec).  For 
oblique  incidence,  the  planar  approximation  is  used  to  reduce  the  problems 
to  two  space  dimensions. 

Results  of  the  initial  portion  of  the  investigation  have  been  previously 
described  in  NRL  Report  6310  (April  1965).  The  present  report  deals 
primarily  v;ith  Impacts  on  relatively  thick  targets ,  and  with  oblique  incidence 
and  hollow  rod  impacts. 

From  the  numerical  solutions,  detailed  plots  of  the  shock  systems  and 
flow  fields  occurring  during  perforation  are  obtained.  From  the  flow  fields 
for  impacts  where  the  ratio  of  target  thickness  to  rod  diameter  is  two  or 
l-eater,  it  is  found  that  the  processes  causing  rod  loss  can  be  divided  into 
three  phases;  (1)  the  Initial  transient  phase  during  which  the  intense  shock 
system  caused  by  impact  is  attenuated  by  rarefactions  from  free  surfaces-, 

(2)  the  quasi-steady  state  phase,  during  which  the  flow  field  remains 
essentially  unchanged  relative  to  a  stagnation  point  near  the  tip  of  the 
advancing  rod  shaft,  and  (3)  the  breakthrough  phase,  during  which  the  rod 
moves  through  the  remaining  shattered  target  material.  Most  of  the  rod  loss 
occurs  during  the  steady  state  phase. 

Quantitative  information  from  the  flow  fields  is  used  to  derive 
relationships  for  predicting  rod  loss  as  a  function  of  velocity,  target  thick¬ 
ness,  projectile  and  target  materials,  and  impact  obliquity.  Comparison 
of  these  relationships  with  experimental  data  obtained  by  NRL  provides 
verification  of  their  validity. 


PROBLEM  STATUS 


This  is  a  final  report. 


AUTHORIZATION 


NRL  Problem  F04-16 
Project  MIPR  PG-3-15 

Office  of  Naval  Research  Contract  NONR  4221(00)(X) 
Air  Force  Project  9850 


Report  prepared  by  Shock  Hydrodynamics,  Inc., 
for  the  Naval  Research  Laboratory,  the  project  being 
sponsored  by  the  Air  Force  Armament  Laboratory, 
AFATL,  Eglin  AFB,  Florida. 


Final  composition,  printing,  and  distribution  of 
this  report  have  been  carried  out  by  the  Naval  Re¬ 
search  Laboratory. 


Manuscript  submitted  May  6,  1968, 


NOTATION 


A,  B  constants  in  £q.  (2-21)  and  Table  2-1 
C  local  sound  speed 

D  rod  diameter  (where  projectile  is  a  rod)  or  plete 
thickness  (where  projectile  is  a  plate) 

e  specific  internal  energy 

e-ji  specific  internal  energy  in  target  material 

ep  specific  internal  energy  in  projectile  material 

ID  inside  diameter  of  hollow  rod 

AL  rod  length  lost  during  perforation 

OD  outside  diameter  of  hollow  rod 

P  pressure 

S  shock  speed  in  undisturbed  material  corresponding 

to  impact  velocity  V 

t  time 

t{  time  required  for  Yp  to  equal  T 

T  target  thickness 

u  particle  velocity 

V  impact  velocity 

Voo  stagnation  point  velocity 

Yp  penetration  (or  open  hole)  depth 

Yg  depth  of  stagnation  point 

e  angle  of  impact  incidence  (measured  between 

impact  velocity  and  surface  normal) 

p  density 

p-j.  target  density 

pp  projectile  density 
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NOTATION  (Continued) 


T  time  characteristic  of  initial  transient  phase 
(in  axisymmetric  impact) 

Tq  time  at  which  initial  one-dimensional  nature  of 
shock  system  is  completely  destroyed  by  lateral 
rarefactions 

Tpj^  time  characteristic  of  Initial  transient  phase  in 
planar  impact 

time  characteristic  of  initial  transient  phase  in 
oblique  Impact  of  rod 


FINAL  REPORT 


HYPERVELOCITY  IMPACT  OF  RODS  ON  FINITE  TARGETS 


1.  INTRODUCTION 

This  is  the  final  report  of  an  analytical  research  investigation  of 
impacts  of  end-oriented  rods  into  finite  targets.  The  study,  based  upon 
two-dimensional  hydrodynamic  PIC  {Particle -in-C ell)  code  solutions  of 
selected  impact  cases,  was  aimed  at  providing  (a)  increased  understanding 
of  the  physical  processes  operating  in  such  impacts,  and  (b)  techniques  for 
predicting  impact  effects  for  different  materials  under  various  conditions  of 
impacting  velocities  and  obliquities.  The  research  was  started  in  19  64.  A 
related  experimental  Investigation  has  been  conducted  during  the  same 
period  by  the  Naval  Research  Laboratory. 

In  the  initial  phase  of  the  study,  emphasis  was  given  to  impacts 
against  relatively  thin  targets.  A  previously-issued  report^  describes 
results  of  the  initial  phase,  and  only  those  data  which  are  necessary  to 
provide  continuity  will  be  repeated  in  this  final  report. 

In  the  final  phase,  which  is  described,  in  detail  in  this  report,  thicker 
targets  were  considered,  along  with  hollow  rods  and  also  oblique  impacts. 


1.1  APPROACH 

The  approach  taken  was  to  use  hydrodynamic  solutions  to  obtain 
detailed  data  regarding  mass  positions,  pressures,  and  particle  velocities 
as  functions  of  time.  These  data  were  analyzed  and  compared,  and  from 
them,  relationships  for  rod  loss  as  a  function  of  Impact  velocity  and  pro¬ 
jectile  target  materials  were  developed. 


1.2  IMPACT  CASES  CONSIDERED 

In  the  overall  program,  sixteen  specific  impact  cases  were  examined 
in  detail.  These  are  summarized  in  Figures  1-1  and  1-2.  Eight  of  these 
were  considered  in  the  initial  phase  of  the  study,  and  are  reported  upon  in 
reference  1. 

Solution  plots  for  some  of  the  problems  were  prepared  in  the  form  of 
movies  which  have  been  delivered  to  the  Naval  Research  Laboratory. 
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jure  1-1  Conditions  For  General  Solid  Rod  Impact  Cases 


Conditions 
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1.3  HYDRODYNAMIC  MODEL 

All  of  the  analyses  performed  during  this  study  assumed  that  the  pro¬ 
jectile  and  target  material  properties  could  be  described  using  a  hydrody¬ 
namic  model.  For  the  velocities  which  were  considered,  this  assumption  is 
valid  for  describing  the  early-stage  effects  and  subsequent  steady  state 
phenomena  in  thicker  plate  impacts.  Later  response,  including  establish¬ 
ment  of  hole  size  and  the  rod  loss,  clearly  involves  material  strength  effects. 
At  the  time  the  study  was  undertaken,  analytical  techniques  incorporating 
elastic -plastic  properties  were  not  yet  available.  It  was  therefore  necessary 
to  extract  criteria  for  rod  loss  from  comparison  with  experiments  and  by 
simplified  analytical  means. 

It  is  clear,  however,  that  an  elastic-plastic  model  is  needed  to 
rigorously  describe  the  final  states  of  response.  Under  other  contracts  , 
comprehensive  hydrodynamic-elastic-plastic  numerical  methods  have  been 
developed  and  future  studies  of  rod  impact  could  employ  such  techniques.  ^  ^ 
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2.  SOLID  ROD  IMPACTS 

In  this  section,  mechanisms  causing  projectile  rod  loss  in  the  normal 
impact  of  solid  rods  are  examined  and  illustrated  by  means  of  plots  from  the 
numerical  solutions.  A  relationship  is  then  derived  for  predicting  rod  loss 
for  various  impacts.  This  relationship  is  compared  with  both  numerical  and 
experimental  results. 


2.1  MECHANISMS  IN  SOLID  ROD  IMPACTS  ON  FINITE  TARGETS 

In  high  velocity  impacts  ,  extreme  pressures  are  created  in  the  pro¬ 
jectile  and  in  the  impacted  region  of  the  target,  and  these  high  pressure 
regimes  spread  as  shocks  which  propagate  outward  to  engulf  adjacent 
material.  The  response  of  the  projectile  and  target  materials  to  this  high 
dynamic  pressure  field  can  be  divided  into  three  phases,  which  we  will 
refer  to  as  the  Initial  transient  phase,  steady  state  phase,  and  breakthrough 
phase. 


2.1.1  Initial  Transient  Phase 

Immediately  upon  impact,  shockwaves  propagate  from  the  interface 
into  both  the  projectile  and  the  target.  Where  the  projectile  face  is  flat, 
the  region  between  the  diverging  shocks  is,  for  a  short  time,  at  a  constant, 
high  pressure  state,  in  which  all  motions  are  one-dimensional  (parallel  with 
the  Impact  velocity).  This  constant  one-dimensional  state  persists  until 
rarefactions  from  free  surfaces  arrive  to  relieve  the  pressure.  Sources  for 
such  rarefactions  are  the  front  and  back  surfaces  of  the  target  plate,  and 
the  sides  of  the  projectile.  (For  long  rod  projectiles,  the  trailing  surface 
of  the  rod  is  too  remote  to  be  a  source  of  rarefactions  during  the  initial 
transient  phase.) 

Where  the  target  is  thin  relative  to  the  projectile  diameter,  the  back 
surface  of  the  target  will  be  directly  accelerated  by  the  one -dimensional 
shock,  and  one-dimensional  pressure  relief,  which  we  term  the  frontal 
rarefaction,  will  propagate  back  into  the  target.  The  effect  of  this  frontal 
rarefaction  on  rod  loss  is  important  in  impacts  on  very  thin  targets,  as  is 
discussed  in  detail  in  reference  1.  As  target  thickness  increases,  however, 
lateral  rarefactions  emanating  from  the  sides  of  the  rod  projectile  and  from 
the  uns hocked  adjacent  target  material  become  the  initial  and  primary  mode 
of  pressure  relief  in  the  shocked  region.  These  rarefactions  converge 
towards  the  axis  of  the  problem ,  and  also  overtake  the  shocks  moving  into 
the  target  and  projectile. 

Figure  2-1  Illustrates  these  shock  systems  using  spatial  plots 
obtained  from  numerical  solutions  in  this  program.  These  plots  show 
pressure  contours  (isobars)  at  selected  intervals.  The  left  boundary  is  an 
axis  of  symmetry,  and  corresponds  to  the  center  line  of  the  rod  projectile. 


(uio)  q: 


Ml 

innnrwinnuHHnMuni  . 

.  jimaaimnaiiBnniUBnmuMHBi 

gg55ggSag58g5Sa5gg^gg8^ 

■■■■■KnivaBaRaaBiiissnniBBraaaiii 

zsmmmmwMK^smmmKmwiammmmmmmmnt 

IBflBB!?fiBBBBBBDBBKA'%9nB0BBBIll 

_ JBBBBBBBnBRBBBBBBIKJIfiBISIISSBBBI 

SSaESSSSSBSaa9SS=B599l5eSIBBBBKBB&^3SSBIS!6BBI 
SSBSBBSSSBSSBBBSSSSBSBBBilfiBaKaiRnKBBBBII  BBI 

BBBBBBBBBBBBBBBBBBIflBBBBBBBBBBBKIBllBBBBBBI 

aaaaaaaaaaaaaaaaeaaaBSeafiBBBBBBissBBiinBBBi 

■■■■■■■■■■■■BQSaaaBilBSBRBilBR9BBB!l£Bi^.l!fBBBI 

BaaaaaaaaaaaafiaBBaaaaBtiBBBBBflBSBKis^^jBBBBi 

5SBBSBBasaSB.iiiSiBBS5a:<Ba!ll!9BBBBBBBBaEi!«BBBBI 

B~BBBBMBBBfiafiBBBBSBfiaBBBflBKilBaBBBR'^^:aBBBBI 
BBBBBBBBBBBBBBBBBBBKBUBBBBBBi!ae:^iiigaiBBBBI 
aaaBaaaaBBaaaBBeBaaaa^asaBRassssBRBBBBBBi 


IBUBBBUBBBBUBBBBUBBBBUBBBBBBBBBBBBBBBBBBBI 


Radius  (cm) 


Biimimimiiiiiiii  _ _ _ _ _ 

aasgggSgagassaSBBBMBHiaiililiaaiBBBBBl 
^gaaliaaBgieaiBgiBaasgggggggggggggggi 
■KfifiBgallBiilBBiOggBggHgHBBggl 
VTBBBSBBBISBinrjIBBBI 


Head  of  Lateral 
Rarefaction  Fan 
At  26,2  ^isec 


aBBBBBEliBBBBBBBBBBBWSBBBBBBBBBBlIBB 

BBBBB!!SiriBflBBRBflBIIR3IRifSllBBSfilBBflBBBIIBB 
■BIBBiiBliBniifilKSIBBWISSinBBSSBklBaiBBBBBlIBB 
BBBBflViBBIllBB9IBBflBIB9BBB8i&Sair#«a!9igBBBf 
BBBBailBBBBKBBBBBBnBBBBRSBSiBBSIBBBBI 

IBKffiSaBBBKgBBSBBBKBlIBBBI 
IBIIIIinBBSIBSBBIiniSSBBlIBBBI 
NIMBBBBBr 


:  Frontal  Rarefaction 


I  I  I  I  I  f 


iBBSsSaEIBBBI 


27 

Radius  (cm) 


Figure  2-1  Development  of  Shock  and  Rarefaction 
Systems  In  Thin  Plate  Impacts 


(A  more  thorough  description  of  the  numerical  technique  and  the  plotting 
formats  is  given  in  Appendix  A.) 

In  Figure  2- la,  diverging  shocks  are  propagating  up  into  the  pro¬ 
jectile  and  also  down  into  the  target.  The  head  of  a  rarefaction  fan  is 
propagating  in  from  the  side.  At  the  somewhat  later  time  represented  in 
Figure  2-lb,  the  shock  propagating  into  the  projectile  remains,  but  the 
shock  in  the  target  has  already  reached  the  rear  surface,  and  a  frontal 
rarefaction  is  propagating  back  into  the  target.  The  lateral  rarefaction  is 
converging  into  the  remaining  high  pressure  region,  and  will  soon  thereafter 
reach  the  axis.  Subsequently,  the  lateral  rarefaction  will  also  overtake  the 
upper  shock  in  the  target,  thus  terminating  the  initial  one-dimensional  flow. 

In  the  example  shown  in  Figure  2-1,  the  one-dimensional  shock  in 
the  target  reached  the  rear  surface  prior  to  being  overtaken  by  the  lateral 
rarefaction.  The  resulting  frontal  rarefaction  was  therefore  also  one- 
dimensional.  This  will  always  be  the  case  for  very  thin  plate  impacts. 

(i.e. ,  T/D  <~  1/2  for  Al-Al  impacts  at  5-10  km/sec.)  For  thicker  targets , 
however,  lateral  rarefactions  will  overtake  the  target  shock  before  it  reaches 
the  rear  surface. 

Termination  of  the  one-dimensional  high  pressure  flow  by  encroach¬ 
ment  of  the  head  of  the  lateral  rarefactions  constitutes,  of  course,  only  the 
start  of  the  pressure  relaxation  process  which  continues  in  the  rarefaction 
fan  until  equilibrium  is  reached.  For  thin  targets,  the  equilibrium  state  is 
at  zero  pressure,  and  this  is  reached  directly  from  the  actions  of  the  frontal 
and  lateral  rarefactions.  In  thicker  targets,  the  prs-ssure  cannot  drop  to 
zero,  but  rather  it  drops  to  a  condition  established  by  material  flow  near 
the  tip  of  the  projectile  as  it  continues  to  advance  into  the  target. 

In  summary,  an  initial  transient  stage  occurs  upon  Impact  during 
which  a  very  intense  shock  system  is  established.  Free  surface  rarefac¬ 
tions  subsequently  attenuate  the  shock  system  and  return  the  system  to 
equilibrium,  thereby  completing  the  transient  phase.  In  thin  targets,  this 
transient  stage  encompasses  virtually  all  of  the  significant  response 
phenomena,  and  it  therefore  completes  the  Impact  process.  In  thick  targets, 
the  transient  stage  leads  into  the  quasi-steady  state  phase  described  in 
the  following  section. 


2.1.2  Steady  State  Phase 

In  thick  finite  targets  (or  in  seml-lnflnlte  targets)  the  transient  phase 
gives  way  to  a  steady  state  phase  during  which  a  deep  hole  is  being  formed 
in  the  target  by  the  advancing  tip  of  the  relatively  narrow  rod.  The  rod  tip 
is  continually  receding  (or  being  worn  away)  during  this  process,  leading 
to  rod  loss.  Impacts  into  thicker  targets  constituted  a  major  area  of  study 
during  the  portion  of  the  program  covered  by  this  report.  Hence,  the 
Identification  and  understanding  of  the  steady  state  processes  are 
emphasized. 
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Figure  2-2  illustrates  some  aspects  of  the  steady  state  phase.  (Again, 
this  is  a  mass  position  plot  from  one  of  the  numerical  solutions.)  The 
undamaged  shaft  of  the  rod  continues  to  travel  into  the  plate  at  essentially 
the  impact  velocity  (4.6  km/sec  in  this  case)  throughout  the  penetration 
process.  At  the  rod-target  Interface,  material  is  being  diverted,  first  in 
the  radial  direction,  and  then  upward,  resulting  in  the  formation  of  an  open 
hole  or  crater.  The  rear  surface  of  the  target  plate  has  been  accelerated  by 
arrival  of  the  initial  impact  shock,  and  a  bulge  is  forming  on  that  surface. 
The  velocity  imparted  to  this  material,  however,  is  smaller  than  the  Impact 
velocity.  The  advancing  rod  is  therefore  impinging  into  slower  moving 
target  material. 

Note  in  Figure  2-2  that  we  define  the  depth  of  the  open  hole  below  the 
original  target  surface  as  the  penetration  depth  Yp.  Evidence  that  an 
essentially  steady  state  condition  is  reached  during  penetration  of  a  thick, 
albeit  finite  target  is  seen  in  the  plot  in  Figure  2-3  of  Yp  vs  time.  The 
values  of  Yp  in  this  plot  were  obtained  directly  from  the  numerical  solution 
of  the  T/D  =  3,  Al-Al  impact  at  4.6  km/sec.  Following  the  initial  transient 
stage,  it  is  seen  that  the  rate  at  which  the  depth  increases,  dYp/dt,  remains 
constant. 

In  steady  state  flow,  we  would  also  expect  to  find  a  stagnation  point 
on  the  axis  just  below  the  bottom  of  the  open  hole.  We  will  define  the  depth 
of  this  point  below  the  original  surface  as  Ys.  The  velocity  of  the  stagna¬ 
tion  point  can  be  evaluated  by  Bernouli's  equation  for  compressible  flow, 

l.e. 


P  1  2 

— +•5-  u  +  e  =  constant  on  a  steady  state  streamline 
^  (particle  path)  (2-1) 

where 

P  =  pressure 
p  =  density 
u  =  particle  velocity 
e  =  specific  internal  energy 

In  a  frame  of  reference  which  moves  at  the  stagnation  point  velocity, 
Vm,  consider  the  motions  of  as -yet  undisturbed  particles  A  and  B  as  they 
approach  the  stagnation  point. 


Depth  (cm) 
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2-2  Definition  of  Penetration  Crater  and  Penetration  Depth 


A 


undisturbed  particle  in 
projectile 
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undisturbed  particle 
in  target 


At  point  A,  P,  ep  =  (land  u  =  V  -  Vc.  Hence,  eqn.  (2-1)  has  a  con¬ 
stant  value  of  1/2  (V  -  on  a  steady  state  streamline  in  the  projectile. 
(V  is  the  Impact  velocity,  and  the  subscripts  P  and  T  refer  to  projectile  and 
target  conditions.) 

At  point  B,  P,ex  =  0  and  u  =  -V«,.  On  a  streamline  in  the  target, 
eqn.  (2-1)  therefore  has  a  constant  value  of  1/2 

Evaluation  of  eqn.  (2-1)  at  the  stagnation  point  leads  to; 

(in  projectile) 

I-  +  ep  =  i  (V  -  V,)^  <2-2) 

Pp  2 


(in  target) 


2 


(2-3) 


Eliminating  P  gives; 

(V  -  Voo)^  Pp  -  P^  =  2(eppp  -  e^p^;,)  (2-4; 

For  similar  material  impacts ,  pp  =  p.j.  and  ep  =  .  In  such  impacts  ,  we 

thus  obtain 


V 


(2-5) 
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Under  conditions  approaching  incompressible  flow,  the  right  hand  side 
of  eqn.  (2-4)  becomes  small  compared  to  either  of  the  terms  on  the  left  side, 
and  vanishes  in  the  limit.  This  permits  us  to  evaluate  the  stagnation  point 
velocity  for  dissimilar  material  impacts,  i.e. 

(V  -  V„)2  Pp  =  m}  Pj 


(2-6) 
Pp 

In  Figure  2-3,  we  have  also  plotted  the  depth  of  the  stagnation  point, 
Yg ,  vs  time.  The  position  of  the  stagnation  point  is  not  readily  apparent 
from  plots  of  numerical  solutions  where  the  laboratory  frame  of  reference  is 
used.  To  locate  the  point,  we  therefore  used  the  numerical  solutions  to 
determine  velocity  vs  depth  at  various  times  along  the  axis ,  and  find  the 
depth  at  which  the  velocity  drops  to  half  the  impact  velocity,  i.e. 


1+y^ 


Velocity 

Along 

Axis 


The  Yg  vs  t  points  in  Figure  2-3  can  be  reasonably  well  fit  by  a  straight 

line,  the  slope  of  which  is  2.3,  or  V/2.  Such  a  fit  is  parallel  to  the  Yp  vst 
data,  and  is  offset  by  about  D/2.  Thus  Yg  =  Yp  +  D/2. 

The  conclusion  of  the  steady  state  phase  is  brought  about  by  the 
interactions  with  the  rear  free  surface.  The  presence  of  this  free  surface 
is  felt,  of  course,  at  an  early  time  in  the  solution,  as  soon  as  the  initial 
shock  first  reaches  that  surface.  To  the  extent  that  this  shock  is  one¬ 
dimensional,  the  free  surface  immediately  achieves  a  velocity  nearly  equal 
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to  the  impact  velocity.  For  the  thicker  targets  which  we  have  been  con¬ 
sidering,  however,  lateral  rarefactions  so  reduce  the  strength  of  the  impulse 
delivered  to  the  rear  surface  that  its  initial  velocity  is  a  fraction  of  the 
impact  velocity.  For  the  T/D  =  2  and  T/D  =  3  targets,  for  example,  the 
initial  velocity  of  the  rear  surface  is  less  than  half  the  impact  velocity. 

The  intermediate  material  between  the  rear  surface  and  the  stagnation  point 
is  at  an  even  lower  velocity.'  As  the  stagnation  point  moves  further  into 
the  target  and  becomes  correspondingly  closer  to  the  rear  free  surface,  the 
elevated  pressures  will  provide  additional  accelerating  force  to  the  inter¬ 
vening  material.  Eventually  all  that  material  will  be  moving  at  the  stagna¬ 
tion  point  velocity  or  higher,  a  condition  which  we  find  to  exist  by  the  time 
Yp  =  T.  The  steady  state  condition  is  assumed  to  be  concluded  at  that  point. 


2.1.3  Breakthrough  Phase 

There  follows  the  final,  or  breakthrough  phase  of  the  impact  response, 
during  which  the  rod  interaction  with  the  target  is  complete  and  the  residual 
length  passes  through  an  open  hole.  It  is  not  practical  to  analyze  this  phase 
in  detail  by  means  of  numerical  hydrodynamic  solutions.  We  believe,  how¬ 
ever,  that  breakthrough  constitutes  a  relatively  unimportant  phase  in  the 
overall  rod  loss  process. 

During  the  breakthrough  phase,  the  rod  must  clear  its  way  through  the 
remaining  rear  surface  debris.  For  Yp  T,  it  is  seen  in  the  numerical  solu¬ 
tions  that  the  average  density  of  such  debris  is  low,  indicating  a  shattered 
state.  The  total  mass  ahead  of  the  projectile  is  therefore  small.  Further¬ 
more,  although  the  magnitude  of  the  debris  velocity  is  less  than  that  of  the 
residual  rods,  some  of  the  debris  is  diverging,  so  that  much  of  the  shattered 
mass  will  move  radially  out  of  the  path  of  the  residual  rod.  We  conclude 
that  the  extent  of  rod  loss  caused  by  further  interactions  between  debris  and 
rod  will  be  small  by  comparison  with  the  rod  loss  due  to  the  initial  transient 
and  to  steady  state  flow.  ^ 


2.2  DEVELOPMENT  QF-*B';QUATI0NS  FOR  PREDICTING  ROD  LOSS 

Further  analysis  of  the  above-described  mechanisms  of  the  transient 
and  steady  state  portions  of  the  penetration  process  permits  the  derivation 
of  a  relationship  for  rod  loss  in  terms  of  impact  velocity,  rod  diameter  and 
material,  and  plate  thickness  and  material.  In  this  derivation,  it  is 
hypothesized  that  no  significant  rod  loss  occurs  after  the  open  hole  depth 
reaches  the  original  rear  surface  of  the  target,  i.e.  ,  when  'Yp  =  T.  The 

validity  of  this  hypothesis  is  examined  further  in  Section  2.3.2,  wherein 
the  rod  loss  equation  is  compared  with  experimental  observations. 
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2,2.1  Basic  Rod  Loss  Equation 

As  noted  in  Section  2.1.2,  at  the  stagnation  point  near  the  bottom  of  the 
hole,  rod  material  is  diverted  radially.  Pressures  above  the  stagnation 
point  drop  rapidly,  so  it  is  reasonable  to  consider  that  the  residual  tip  of 
the  intact  rod  is  always  located  at  the  bottom  of  the  hole.  Thus,  if  rod  loss 
is  completed  when  Yp  =  T,  the  final  .ip  of  the  rod  will  be  located  at  the 
position  of  the  original  rear  surface  of  the  plate,  as  indicated  in  Figure  2-4b. 
At  the  time  of  impac  chis  final  rod  tip  is  located,  as  shown  in  Figure  2-4a, 
at  a  point  AL  ahead  of  the  target  front  surface,  where  AL  is  the  length  of 
rod  loss. 

Defining  tj  as  the  time  required  for  the  open  hole  depth,  Y„,  to  equal 
the  plate  thickness,  T,  it  is  seen  that 

Vtf=  AL+  T  (2-7) 

Now  consider  the  manner  in  which  the  penetration  depth  Yp  varies 
with  time.  As  was  seen  in  Section  2,1.2,  the  steady  state  penetration  rate 
for  similar  material  impacts  is  V/2.  In  Figure  2-3,  the  plot  of  Yp  vs  t 
approaches  V/2  asymtotically.  This  suggests  an  exponential  fit  of  the  form 


Yp  =  y  [t  -  T  d-e"^/'^)  }  (2-8) 

in  which  T  is  the  time  characteristic  of  the  transient  phase.  For  the 
general  case  of  non-similar  Impacts,  V®  would  be  substituted  for  V/2, 
giving 

Yp  =  V®  {t  -  T  (1-e”^^'^)  }  (2-9) 

Since  we  assume  that  rod  loss  stops  when  Yp  =  T  at  time  t£,  eqn. 
(2-8)  becomes 


T  =  I  {tf  -  T  (1-e  } 


or  alternatively 


I^VT  ril  _  l  + 

D  2d'-t  ^ 


Substituting  t^  from  eqn,  (2-7), 

Al=vt^t  _;^  -(^^) 
D  D  D  D 


(2-10) 


(2-11) 


(2-12) 


view  at  t 
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Our  assumption  that  rod  loss  stops  at  Yp  =  T  is  limited  to  large  values  of 
T/D.  For  these  cases,  the  exponential  term  in  eqn.  (2-12)  tecomes  insigni¬ 
ficant,  and  we  obtain  the  basic  relationship  for  rod  loss  in  large  T/D  Impact 


cases . 


vr  +1 

D  D  D 


(2-13) 


2.2.2  Evaluation  of  Vt/D  From  Numerical  Data 

As  t  -*  ®  ,  the  exponential  relationship  for  Yp  vs  t  in  eqn.  (2-8) 
becomes 


Yp=  ^  (t  -  T) 


(2-14) 


which  is  the  equation  of  the  straight  line  to  which  eqn.  (2-  *s  asymptotic. 
Equations  of  this  form  can  be  nominally  fit  to  Yp  vs  t  data  obtained  from 
numerical  solutions  to  determine  approximate  values  for  time  characteristics. 
This  is  done  in  Figure  2-5  for  im.pact  Case  11.  The  time  intercept  of  the 
V/2  asymptote  gives  T  to  be  200.  Evidence  that  both  this  value  and  the 
hypothesized  exponential  relationship  in  eqn.  (2-8)  are  valid  is  shown  by 
the  excellent  fit  of  the  numerical  points  in  Figure  2-5. 

For  aluminum-aluminum  impacts  at  4.6  km/sec,  the  value  of  Vt/D  in 
eqn.  (2-13)  is  thus  established  as  .5,  and  that  relationship  becomes 


K 

D  D 


(2-15) 


Note  that  although  D  =  60  cm  was  used  in  the  numerical  solutions  (and  in 
evaluating  t)  Vt/D  Is  dimensionless  and  can  be  applied  to  any  scale 
impact. 


2.2.3  Vt/D  For  Other  Similar-Material  Impacts  and  Velocities 

T  is  related  to  the  time  required  for  the  one-dimensional  nature  of  the 
initial  shock  systems  to  be  destroyed  by  the  lateral  rarefaction  waves, 
which  we  will  call  Tq.  The  relationship  between  T  and  Tq  is  probably 
complex,  but  for  the  present  purposes  it  is  reasonable  to  assume  that  t  is 
proportional  to  T^,  i.e. 


T  =  aT, 


(2-16) 


From  this ,  it  is  possible  to  establish  the  value  of  Vt/D  for  other  conditions 
of  impacting  materials  and  velocities. 


J  'It  4 
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The  head  of  the  lateral  rarefaction  fan  propagates  Into  the  initial  one- 
dimensional  high  pressure  region  with  local  sound  speed  C.  Thus 


’’o  = 


D/2 


J -  (s -y/iy 


(2-17) 


where  S  Is  the  shock  speed  corresponding  to  impact  velocity  V.  C  and  S 
can  be  evaluated  from  an  equation  of  state.  Using  the  Los  Alamos  equation 
of  state,  Tq  for  aluminum  Impacts  at  4.6  km/nsec  is  54.25,  This  leads  to 
a  proportionality  constant  in  eqn.  (2-16)  of  3.69.  Thus  we  arrive  at  the 
following  general  relationship  for  evaluating  Vt/D  in  eqn,  (2-13)  for  similar 
material  Impacts: 


W  ^  1,84  V 

D  yc2-  (S-V/2)^ 


(2-18) 


2.2.4  Simple  Method  for  Estimating  Vt/D  Term 

Where  it  is  inconvenient  to  obtain  the  required  one-dimensional  shock 
Information  required  for  eqn.  (2-18),  (i.e.  C  and  S  as  functions  of  V) ,  a 
less  accurate  but  simpler  method  can  be  used  to  obtain  rod  loss  estimates 
for  a  given  impact.  This  method  is  based  upon  the  two  assumptions 

a)  Shock  velocity  is  linear  with  particle  velocity,  i.e. 

S  =  A  +  Bu  (2-19) 

b)  Sound  velocity  is  related  to  shock  velocity  by 

C'*'  =  .49  S^  +  (S  -  V/2)^  (2-20) 

The  relationship  in  eqn.  (2-20)  is  suggested  in  reference  4. 

Substituting  eqn.  (2-19)  and  (2-20)  into  (2-18),  leads  to  the 
approximation 

Vt  _  2.63  V 

D  'TTW' 

which  can  be  used  in  eqn.  (2-13).  Values  of  A  and  B  in  this  approxima¬ 
tion  are  given  in  Table  2-1,  along  with  calculated  values  of  Vt/D  for 
Impacts  at  4.  6  and  6.  1  km/sec. 


Table  2-1 

Values  of  A,  B,  and  Vt/D  For 
Similar-Material  Impacts  of  Various  Materials 


19 


Material 

Density 

gm/cmi^ 

A 

km/sec 

B 

Reference  Vr/D 

V  =4. 6  km/sec  6 

.  1  km/sec 

Ag 

10.49 

3.215 

1.643 

5 

1.7 

2.0 

A1 

2.785 

5.37 

1.339 

6 

1.4 

1.7 

Au 

19.24 

3.059 

1.608 

5 

1.8 

2.0 

Be 

1.845 

7.975 

1.091 

5 

1.2 

1.4 

Cd 

8.64 

2.408 

1.718 

5 

1.9 

2.1 

Co 

8.82 

4.652 

1.506 

5 

1.5 

1.7 

Cr 

7.13 

5. 176 

1.537 

5 

1.4 

1.6 

Cu 

8.9 

3.972 

1.478 

5 

1.6 

1.9 

Fe 

7.85 

3.8 

1.580 

6 

1.6 

1.9 

In 

7.27 

2. 370 

1.608 

5 

2.0 

2.2 

Mg 

1.735 

4.49  3 

1.266 

5 

1.6 

1.9 

Mo 

10.2 

5.173 

1.204 

5 

1.5 

1.8 

Nb 

8.604 

4.447 

1.212 

5 

1.7 

2.0 

Ni 

8.86 

4.667 

1.410 

5 

1.5 

1.8 

Pb 

11.34 

2.066 

1.517 

5 

2.2 

2.4 

Pd 

11.95 

3.793 

1.922 

5 

1.5 

1.7 

Pt 

21.37 

3.671 

1.405 

5 

1.8 

2.0 

Rh 

12.42 

4.68 

1.645 

5 

1.4 

1.7 

Sn 

7.28 

2.668 

1.428 

5 

2.0 

2.3 

Ta 

16.46 

3.  374 

1.155 

5 

2.0 

2.3 

Th 

11.  68 

2.079 

1.381 

5 

2.3 

2.6 

Ti 

4.51 

4.786 

1.066 

5 

1.7 

2.0 

T1 

11.84 

1.821 

1.566 

5 

2.2 

2.4 

W 

19.  17 

4.005 

1.268 

6 

1.7 

2.0 

Zn 

7. 135 

3.042 

1.576 

5 

1.8 

2.0 

Zr 

6.49 

3.771 

.933 

5 

2.0 

2.4 
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Indications  of  the  accuracy  of  this  approximation  are  seen  in  Figures 
2-6  and  2-7  for  A1  vs  A1  and  Fe  vs  Fe  impacts.  These  plots  compare  Vt/D 
as  calculated  by  the  more  exact  eqn.  (2-18)  with  the  approximation  of  eqn. 
(2-21).  It  is  seen  that  eqn.  (2-21)  understates  Vr/D  for  aluminum  by 
approximately  0.1  -  0.2  diameters  over  a  wide  velocity  range.  In  iron, 
eqn.  (2-21)  overstates  Vt/D  by  up  to  0.3  diameters  at  velocities  above 
0.2  cm/i-isec.  While  these  discrepancies  are  relatively  small,  some  changes 
in  the  ordering  of  material  effectiveness  may  occur  when  the  different  rod 
loss  prediction  equations  are  used. 


2.2.5  Extension  to  Dissimilar-Material  Impacts 

An  insufficient  number  of  dissimilar-material  cases  were  solved 
numerically  to  permit  us  to  directly  derive  relationships  for  rod  loss  in  such 
cases.  However,  we  believe  that  adequate  predictions  can  be  made  for 
dissimilar  Impacts  by  using  incompressible  jet  theory  to  modify  eqn.  (2-13). 
Thus 


-  T 


At 

D 


) 


T  _  X 


(2-22) 


First  calculate  Al/D  from  eqn.  (2-13)  for  an  impact  where  the  projectile  is 
assumed  to  be  of  the  same  material  as  the  actual  target  material.  Then 
multiply  this  result  by  the  square  root  of  the  ratio  of  the  densities  of  the 
target  material  and  the  actual  projectile  material  to  obtain  the  predicted 
rod  loss. 


2.3  COMPARISON  OF  ROD  LOSS  EQUATION  WITH  NUMERICAL 
SOLUTIONS  AND  EXPERIMENTAL  DATA 


2.3,1  Numerical  Solutions 

The  detailed  results  of  the  numerical  solutions  conducted  during  this 
investigation  are  given  by  spatial  plots  of  mass  positions,  velocity  fields, 
and  pressure  contours  in  reference  1  and  in  the  Appendix  to  this  report. 

From  these  solutions,  penetration  depths  (or  the  depth  of  the  open  hole, 

Yp)  can  be  plotted  as  functions  of  time.  To  obtain  the  penetration  depth  from 
the  numerical  solutions,  density  contours  near  the  bottom  of  the  open  hole 
were  examined.  The  density  in  the  hole  is  of  course  zero,  but  the  walls 
and  bottoms  of  the  holes  are  seen  as  sharp  density  gradients  in  the  numeri¬ 
cal  results,  rather  than  as  discontinuities.  For  consistency,  we  define 
the  surface  as  the  density  contour  where  the  compression  is  0.5. 

Plots  of  penetration  depth  vs  time  appear  in  Figure  2-8  for  Al-Al 
impacts  at  4,6  km/sec  into  targets  having  T/D  values  of  1/4,  1/2,  1,2, 

3  and  5.  For  the  small  T/D  Impacts,  the  penetration  rates,  dYp/dt 
accelerate  through  the  steady  state  velocity  Voo  =  V/2  due  to  the  early 
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CONFIDENTIAL 


influence  of  the  rear  free  surface.  For  the  T/D  =  2  and  3  cases  y  the  pene¬ 
tration  rates  asymtotically  approach  V/2 .  {Case  12,  T/D  =  5,  produced 
somewhat  erratic  Yp  vs  t  results  due  to  the  numerous  grid  changes  necessi¬ 
tated  by  the  long  running  time.) 


2.3.2  Comparisons  With  Experimental  Data  (Confidential) 

I 

A  series  of  tests  have  been  performed  by  the  Naval  Research  Laboratory 
in  a  parallel  experimental  Investigation  of  rod  impact  phenomena.  Primary 
observations  in  these  experiments  were  the  residual  rod  lengths  after  per¬ 
forating  plates  of  various  thicknesses.  From  these  observations  the  length 
of  rod  shattered  as  a  function  of  plate  thickness  was  determined. 


2 . 3 . 2 . 1  Aluminum -Aluminum  Targets  (Confidential) 

Figure  2-9  summarizes  Al-Al  impact  data  at  4.6  km/sec.  NRL  has  fit 
these  data  by  the  straight  dotted  line.  It  is  seen  that  the  numerically- 
derived  rod  loss  equation  predicts  the  experimental  data  to  within  about  0 . 3 
rod  diameters.  Although  experimental  scatter  is  evident,  the  under¬ 
prediction  appears  systematic.  This  suggests  the  magnitude  of  rod  erosion 
which  may  be  attributed  to  the  breakthrough  phase  (which  was  neglected  in 
the  development  of  eqn.  (2-13)  ).  For  a  T/D  =  3  impact,  each  phase  would 


account  for  the  following  los 

ses; 

Initial  Transient: 

(Vt/D) 

1.5 

Steady  State; 

T/D 

3. 

Breakthrough; 

0.3 

0.3 

Thus  a  third  term  might  appropriately  be  appended  to  eqn.  (2-13)  to  account 
for  the  breakthrough  phase,  i.e. 


^  _  Vt  ^  T  ^ 
D  -"D  ^  d"  + 


0.3 


(2-23) 


Only  very  limited  data  are  available  regarding  the  effects  of  velocity 
upon  rod  loss  for  target  thicknesses  equal  to  or  above  the  T/D  ^  2  condi¬ 
tion.  These  are  shown  in  Figure  2-10  and  compared  with  the  prediction  of 
eqn.  (2-23). 

As  previously  emphasized,  the  rod  loss  equations  have  been  derived 
for  relatively  thick  targets  (i.e.  T/D  ^  2).  When  applied  to  thin  targets, 
these  relationships  tend  to  overpredict  the  rod  loss.  To  illustrate, 

T/D  =  1  data  are  in  .luded  in  Figure  2-10,  and  eqn,  (2-23)  is  shown  by  a 
dotted  curve  for  comparison.  The  experimental  data  all  fall  somewhat 
below  the  theoretical  curve. 
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Figure  2-9  Comparison  of  Rod  Loss  Prediction  Equation  With  Experimental  Data  for  Aluminum 
Alumintun  Impacts  at  4.6  Km/sec 
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Velocity  (km/sec) 

Figure  2-10  Comparison  of  Rod  Loss  Prediction  Equation 
with  Experimental  Data  for  Aluminu'.n- 
Alumijium  Impacts  Showing  K.  of  Velocity 
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2 . 3 . 2 . 2  Dissimilar  Targets  (Confidential) 

The  NRL  experiments  also  Included  a  number  of  impacts  of  other 
dissimilar  materials.  Many  of  these  impacts  were  against  T/D  =  1.0 
targets,  and  are  again  below  the  thicknesses  for  which  eqn.  (2-13)  and 
(2-22)  ware  derived.  Figure  2-11  shows  the  experimental  points  and  the 
prediction  curves  for  the  materials  and  velocities  involved  as  determined 
from  the  rod  loss  eqn.  ^2-23)  as  modified  by  eqn.  (2-22)  for  dissimilar 
materials . 


2.4  THIN  PLATE  IMPACTS  (Confidential) 

The  rod  loss  equation  developed  in  Section  2.2.1  applies  to  impacts 
where  the  target  thickness  is  of  the  order  of  twice  the  rod  diameter  or 
greater  (T/D  ^  2) .  For  thinner  targets  ,  the  penetration  phenomena  never 
approach  the  quasi-steady  state  conditions  which  eqn.  (2-13)  presumes. 

During  the  initial  part  of  this  investigation  reported  in  reference  1, 
a  rod  shattering  criterion  was  hypothesized  for  application  to  thin  plate 
impacts.  This  criterion  used  the  peak  pressure  propagating  up  the  rod 
axis  to  establish  the  length  of  rod  shattered  (i.e.  shattering  would  stop 
when  the  peak  dropped  below  a  critical  level.)  By  preliminary  comparisons 
with  experiments,  this  critical  value  was  set  at  100  kb  for  aluminum  alloys. 
Application  of  this  criterion  to  aluminum  rod  impacts  into  thin  aluminum 
plates  led  to  the  conclusion  that  a  plateau  would  exist  in  rod  loss  as  the 
target  thickness  increased.  Based  on  this  study,  the  relationship  of 
AL/D  to  T/D  indicated  by  the  heavy  line  in  Figure  2-12  was  suggested. 
Experimental  data  from  NRL  appears  to  confirm  the  existence  of  such  a 
plateau,  as  is  seen  in  Figure  2-13  (the  heavy  line  is  a  fit  by  NRl) . 

Use  of  the  peak  pressure  criterion,  however,  also  led  to  the  pre¬ 
diction  in  reference  1  that  rod  loss  would  sharply  increase  with  impact 
velocity.  Comparison  of  the  experimental  data  for  4.6  and  6. 1  km/sec 
Impacts  in  Figure  2-13  shows  this  prediction  to  be  incorrect.  Thus  while 
the  peak  pressure  criterion  indicates  the  existence  of  a  rod  shattering 
plateau  which  is  apparently  verified  by  experimental  data,  this  simple 
criterion  is  inadequate  to  describe  the  several  processes  which  are  involve  ! 
in  thin  plate  rod  loss.  To  correctly  describe  these  processes  almost 
certainly  requires  a  comprehensive  model  which  incorporates  material 
strength  effects  as  well  as  fracture  criteria  and  the  properties  of  the 
materials  after  failure . 

In  the  absence  of  a  rigorous  analytical  treatment  of  the  very  thin 
plate  regime,  it  is  of  interest  to  note  that  the  experimental  data  in  this 
regime  can  be  fit  reasonably  well  by  the  exponential  form  given  in  eqn. 
(2-12).  To  this  equation  can  be  added  the  0.3  offset  term  indicated  by 
the  experimental  comparison  in  Section  2. 3.2. 1. 
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NRL  Rod  Loss  Observations  Suggesting  Plateau  for 
Aluminum  -  Aluminum  Impacts 
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Figure  2-14  shows  a  comparison  of  this  equation  with  experimental 
data  over  the  entire  T/D  range  of  available  data.  We  emphasize  that  while 
there  is  a  physical  basis  for  applying  eqn.  (2-12)  to  thick  targets,  it 
represents  only  a  convenient  empirical  fit  of  experimental  data  for  thin 
targets  (T/D  less  than  1). 
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Comparison  of  Experimental  Rod  Loss  Data  With  Exponential  Equation  for  Aluminum 
Aluminum  Impacts  at  4.6  km/sec 
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3.  OBLIQUE  IMPACTS 


3. 1  planar  approximation 

Three  space  dimensions  are  required  to  properly  define  oblique  Impact 
geometries.  At  the  time  the  solutions  in  this  study  were  performed,  no 
suitable  three-dimensional  numerical  technique  was  available.  Invoking 
certain  approximations,  however,  it  is  possible  to  analyze  oblique  impacts 
in  two  space  dimensions,  and  from  such  solutions  to  obtain  semi-quantitative 
information  about  such  impacts. 

For  these  analyses ,  end-oriented  rods  Impacting  at  oblique  incidence 
are  approximated  by  the  two-dimensional  geometry  of  an  infinitely  wide  plate 
Impacting  on  its  edge  against  an  infinitely  wide  target,  as  Illustrated  in 
Figure  3-1.  A  plane  erected  normal  to  the  target  and  to  the  projectile  plate 
produces  a  two-dimensional  cross  section  of  the  problem.  We  will  refer  to 
the  numerical  analyses  of  such  plane  strain  cases  as  planar  solutions. 

The  assumption  of  plane  strain  to  describe  oblique  incidence  precludes 
any  acticn  in  the  direction  normal  to  the  plane  of  the  solution  (z-direction) . 
Thus,  geometric  divergence  of  the  shockwave  in  the  z-direction  is  ignored. 
Conversely,  free  surfaces  do  not  exist  in  the  z-direction;  hence  convergence 
of  rarefactions  from  that  direction  is  ignored.  Certain  qualitative  observa¬ 
tions  can  therefore  be  made  regarding  differences  between  true  three-dimen¬ 
sional  geometries  and  planar  approximations; 

a.  At  a  given  depth  in  the  target,  the  planar  solutions 
will  predict  a  higher  peak  shock  pressure,  energy, 
and  momentum  density  than  exists  in  the  real  three- 
dimensional  case. 

b.  The  initial  transient  stage  will  be  prolonged  in  the 
planar  solution  (as  compared  to  the  three-dimensional 
case) . 

It  is  important  to  note,  however,  that  the  steady  state  characteristics 
(stagnation  pressure,  velocity,  etc.)  are  the  same  for  both  the  planar  and 
three-dimensional  cases.  Thus,  after  the  initial  transient  has  passed,  the 
rate  of  rod  loss  and  of  hole  growth  will  be  the  same  for  both  cases. 


3.2  NUMERICAL  SOLUTIONS 

Planar  solutions  were  obtained  for  30°  and  60°  incidence  Impacts.  A 
normal  incidence  case  was  also  Included  in  order  to  provide  a  direct  com¬ 
parison  between  an  approximated  (planar)  solution  and  a  correct  axisymme- 
trlc  solution  of  the  same  impact.  From  this  comparison,  a  method  was 
devised  for  converting  certain  data  obtained  from  planar  solutions  to  three- 
dimensional  impact  problems. 


Figure  3-1  AxlsiTnmetrlc  vs  Planar  Solutions 


3.2.1  Planar  vs  Axisymmetric  Solutions 

The  Impact  case  chosen  for  this  comparison  was  for  aluminum  vs 
aluminum  at  4.6  km/sec,  where  the  target  thickness  is  twice  the  projectile 
diameter  or  thickness  (i.e.  ,  T/D  =  2.)  This  impact  geometry  was  analyzed 
as  an  ;  ■cisymmetric  problem  in  Case  4  in  reference  1. 

At  early  stages,  the  plotted  solutions  are  r3arly  identical.  In  Figdre 
3-2,  the  positions  of  the  heads  of  the  rarefaction  fans  are  compared  1-3  -  14 
nsec  after  impact.  In  both  cases,  the  fans  emanate  from  the  original  corner 
made  by  the  side  of  the  projectile  and  the  target  face,  and  the  head  of  the 
wave  propagates  into  the  high  pressure  state  at  the  same  velocity.  Behind 
the  front,  however,  the  relaxation  of  pressure  occurs  more  rapidly  in  the 
axisymmetric  case,  due  to  the  diverging  geometry. 

Figures  3-3  and  3-4  compare  the  mass  and  velocity  plots  at  about 
122  jjsec,  by  which  time  significant  differences  are  appearing.  In  Figure 
3-3,  the  depth  of  the  open  hole  is  seen  to  be  appreciably  greater  in  the 
axisymmetric  case.  Conversely,  in  the  velocity  plot  in  Figure  3-4,  a 
somewhat  larger  region  is  seen  to  be  engulfed  by  the  shock  in  the  planar 
case. 


The  depth  of  the  open  hole  continues  to  be  greater  in  the  axisymmetric 
case  throughout  the  solution.  This  is  seen  in  Figure  3-5,  in  which  Yp  vs  t 
data  are  plotted  for  both  cases .  These  curves  are  seen  to  be  nominally 
parallel  and  separated  by  approximately  75  lisec.  (It  takes  the  planar 
solution  75  lisec  longer  than  the  axisymmetric  solutions  to  open  the  same 
hole  depth.)  This  difference  arises  entirely  in  the  initial  transient  stage, 
which  persists  for  a  longer  time  in  the  planar  case  due  to  the  slower  relaxa¬ 
tion  by  lateral  rarefactions.  Thus,  the  time  charuuteristic  rpL  Initial 

transient  in  the  planar  solution  is  275  p.sec,  as  compared  with  t  =  200  psec 
for  the  axisymmetric  solution,  or 


’’’s-D  “  '^PL 


200. 

275^ 


.73  T 


PL 


(3-1) 


Once  the  steady  state  is  reached,  hole  growth  proceeds  at  the  same  rate  for 
both  cases.  Insofar  as  open  hole  growth  curves  are  concerned,  the  planar 
solutions  can  be  converted  to  an  equivalent  three-dimensional  solution  by 
transposing  the  curve  in  accordance  witheqn.  (3-1). 

The  fact  that  t  is  significantly  larger  in  planar  cases  has  another 
effect  on  these  solutions.  As  is  evident  in  Figure  2-5,  steady  state 
conditions  are  not  achieved  until  two  or  three  time  characteristics,  t  , 
have  elapsed.  For  the  impact  treated  in  Case  31,  this  would  mean  that 
500  -  700  pisec  will  elapse  before  the  steady  state  is  reached.  By  this 
time,  however,  the  open  hole  is  approaching  the  rear  target  surface,  and 
no  steady  state  is  possible.  Tills  situation  is  analogous  to  small  T/D 
cases  in  axisymmetric  geometry.  Thus  the  effective  T/D  for  a  planar  solu¬ 
tion  is  smaller  than  indicated  by  the  dimensions  of  the  problem. 
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Figure  3-3  Comparison  of  Mass  Positions  at  122  nsec  for  Axisymmetric  and  Planar 
Solutions  of  Normcil  Impact  of  Aluminum  on  Aluminum  at  4.6  km/sec 
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3-4  Comparison  of  Velocity  Fields  at  122  psec  for  Axisymmetric  and  Planar 
Solutions  of  Normal  Impact  of  Aluminum  on  Aluminum  at  4.6  km/sec 
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A  second  important  quantitative  comparison  between  axisymmetrlc  and 
planar  representation  is  the  peak  pressure  on  the  axis.  Figure  3-6  shows 
these  data.  Large  differences  are  seen,  with  the  plane  strain  case  retaining 
a  given  pressure  level  to  greater  depths  than  the  axisymmetrlc  case.  Again, 
this  is  due  to  the  absence  of  lateral  rarefactions  in  the  planar  case. 

In  summary,  it  is  seen  that  the  planar  representation  of  three- 
dimensional  geometries  has  two  important  effects  on  penetration  and  pro¬ 
jectile  shattering. 

1.  Because  the  transient  phase  is  prolonged  in  planar 
geometry,  the  length  of  rod  lost  during  this  relatively 
inefficient  period  is  greater  than  in  the  case  in  the 
axial  symmetry.  Hence,  a  planar  analysis  will  tend 
to  overestimate  the  actual  projectile  loss. 

2.  Because  the  peak  pressure  reaching  the  rear  surface 
is  greater  in  the  planar  case  (and  also  because  the 
initial  transient  lasts  longer)  a  planar  solution  for  a 
given  T/D  corresponds  to  a  smaller  T/D  for  an 
axisymmetrlc  solution. 


3.2.2  10°  and  60°  Impacts 

Spatial  plots  and  descriptions  of  the  planar  solutions  of  30°  and  60° 
impacts  are  given  in  the  Appendix.  Figure  3-7,  showing  plots  of  the  perie- 
tration  (open  hole)  depth  Yp  vs  time  obtained  from  the  solutions,  provides 
the  basis  for  modifying  the  rod  loss  eqns.  (2-12)  and  (2-13)  for  oblique 
impacts.  Figure  3-7  compares  Yp  vs  t  for  the  30°  and  60°  cases  with  the 
data  from  the  planar  solution  of  a  normal  Incidence  case.  It  is  seen  that 
the  angle  of  obliquity  does  not  significantly  affect  either  the  duration  of 
the  initial  transient  stage,  or  the  steady  state  perforation  rate. 


3.3  MODIFICATION  OF  ROD  LOSS  EQUATION  FOR  OBLIQUE  IMPACTS 

Equation  (2-13)  gives  the  rod  loss  relationship  for  normal  incidence, 

i.e. 


_  Vt  .  X 

D  "  D  D 


(2-13) 


The  transient  and  steady  state  terms  of  this  equation  can  be  appropriately 
modified  for  oblique  impacts  by  reference  to  the  numerical  solutions  and  to 
the  oblique  geometry. 


i 
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Penetration  Depth  vs  Time  from  Planar  Solutions  of  Normal,  30 
•d  60°  Rod  Impacts  (Cases  31,  32,  and  33)  at  4.6  km/sec 
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3.3.1  Trans li^nt  Term 

As  seen  in  Figure  3-7,  impact  obliquity  does  not  significantly  affect 
the  duration  of  the  transient  stage  in  planar  solutions.  Thus  for  30°,  60°, 
and  normal  incidence,  Tp^  =  275  nsec,  and,  using  eqn.  (3-1),  Tgj-j  =  200 nsec. 

We  conclude  that  the  transient  term  Vt/D  in  eqn.  (2-13)  is  Independent  of 
obliquity  angle,  at  least  for  values  up  to  60°. 


3.3.2  Steady  State  Term 

Figure  3-7  also  shows  th^t  the  steady  state  rate  of  hole  growth  is 
independent  of  obliquity.  Howe^'er,  in  oblique  impacts,  the  effective  thick¬ 
ness  of  the  target  is  increased  as  depicted  in  the  following  sketch. 


This  effective  thickness  should  therefore  be  substituted  in  eqn.  (2-13). 
Thus 


AL  _  IM.  i  tan  fl  +  ^ 
D  cos  0  ^  2  ®  ^  D 


(3-2) 


3.4  EXPERIMENTAL  COMPARISON  (Confidential) 

Only  a  very  small  amount  of  experimental  data  exist  for  30°  and  60° 
obliquity  Impacts  for  cases  where  the  T/D  ratio  is  two  or  greater.  Two  such 
points  are  shown  in  Figure  3-8  ,  along  with  straight  lines  representing 
eqn.  (3-2). 
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Figure  3-8  Comparison  of  Rod  Loss  Equation 
with  Experimental  Data  for  Oblique  Incidence 
of  Aluminum  Impacts  into  Aluminiim  at  4.6 
km/ sec 


4.  HOLLOW  ROD  IMPACTS 


The  same  mechanisms  which  have  been  described  for  solid  rod  impacts 
occur  in  hollow  rod  cases,  and  the  same  three  phases  of  impact  are  present. 
The  hollow  geometry  modifies  the  relative  importance  of  these  mechanism  ' 
and  the  duration  of  the  phases,  and  also  introduces  at  least  one  new 
mechanism.  These  differences  are  outlined  in  the  following  section,  and 
their  effects  are  incorporated  into  the  rod  loss  relationship. 


4.1  COMPARISON  OF  MECHANISMS 

In  the  end-on  impact  of  either  a  solid  or  hollow  rod  into  a  finite  target, 
the  pressures  produced  at  the  solid-solid  interfaces  are  as  predicted  by  the 
usual  Hugoniot  jump  conditions.  Where  the  Impacting  surfaces  are  flat, 
the  initial  shock  systems  in  either  case  are  one-dimensional  until  the  effects 
of  lateral  free  surfaces  are  felt.  With  a  solid  rod,  the  active  free  surfaces 
are  the  adjacent  front  of  the  target,  and  the  sides  of  the  rod  projectile.  With 
hollow  rods,  the  inner  surface  of  the  projectile,  as  well  as  the  target  surface 
in  the  center  of  the  impacting  hollow  rod,  are  also  sources  of  pressure  relief. 
By  comparison  with  solid  rods  of  the  same  diameter  or  mass  per  unit  length, 
the  distances  which  rarefaction  waves  must  travel  to  engulf  the  1-d  region 
are  always  shorter  in  hollow  rods.  As  a  consequence,  the  duration  of  the 
initial  transient  stage  is  of  shorter  duration. 

Initially  the  shock  system  and  mass  flow  resulting  from  the  impact  of 
the  wall  of  a  hollow  rod  are  essentially  identical  to  those  which  are  produced 
by  edge-on  impact  of  a  flat  plate.  This  is  seen  in  Figure  4-1.  The  left  axis 
is  the  axis  of  symmetry  of  the  hollow  rod.  Hence  the  projectile  is  repre¬ 
sented  as  the  cross-section  of  one  side  of  the  wall.  At  this  early  stage,  the 
mass  positions  in  the  impact  region  are  similar  on  both  the  inside  and  out¬ 
side  of  the  wall.  Subsequently,  however,  the  shocks  and  mass  flow  on  the 
inner  side  of  the  rod  wall  converge  towards  the  axis.  The  results  of  this 
convergence  begin  to  become  evident  in  Figure  4-2.  The  converging  ejecta 
is  moving  upward  in  the  hollow  of  the  rod  at  a  considerably  higher  velocity 
than  is  the  case  for  the  diverging  ejecta  outside.  Shock  pressures  along 
the  axis  are  in  excess  of  700  kb  (approximately  as  high  as  the  original 
impact  pressure) . 

This  extreme  pressure  is  relieved  by  rarefactions  from  the  free  surface 
at  the  axis ,  causing  upward  acceleration  of  target  material  along  the  axis , 
This  counterflowing  mass  has  relatively  low  density,  but  its  velocity  is  of 
the  same  magnitude  as  the  impact  speed. 

Once  the  elevated  pressures  reach  the  rear  free  surface,  that  surface 
is  of  course  also  accelerated.  Along  the  axis,  mass  is  thus  flowing  both 
upward  and  downward.  At  early  times,  we  therefore  find  that  densities  are 
low  along  the  axis  . 
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figure  4-1  Hollow  Aluminum  Rod  Impact 
At  4.  6  km/sec  into  Aluminum 
At  t  =  14.0  nsec 


Radius  (cm) 


Depth  (cm) 


Following  the  initial  transient  phase,  the  flow  in  hollow  rod  impacts 
enters  the  quasi-steady  state  phase.  The  steady  state  velocity  Vos  will  be 
the  same  as  in  solid  rod  or  planar  cases,  since  the  steady  state  equation 
along  a  stream  line  remains  the  same.  Instead  of  a  stagnation  point,  (as 
in  the  axisymmetric  case)  or  a  stagnation  line  (as  in  the  planar  case)  ,  there 
is  a  stagnation  "ring"  . 


4.2  MODIFICATION  OF  ROD  LOSS  EQUATION  FOR  HOLLOW  RODS 

Equation  (2-13)  gives  the  basic  rod  loss  relationship  for  solid  rods. 

For  hollow  rods,  this  relationship  continues  to  be  valid,  except  that  the 
time  characteristic  of  the  transient  stage,  t,  will  be  reduced  to  reflect  the 
shorter  duration  of  that  phase. 

Figure  4-3  shows  the  penetration  depth  Yp  vs  time  for  two  hollow  rod 
solutions  (ID/OD  =  1/3  and  2/3)  and  compares  these  with  a  conesponding 
solid  rod  case.  The  shorter  duration  of  the  initial  transient  stage  is  apparent: 
T  is  approximately  125  for  both  of  the  hollow  rod  cases,  as  compared  to  200 
for  the  solid  rod.  The  use  of  T  =  125  in  the  exponential  eqn.  (2-8)  for  Yp  vs 
t  produces  the  solid  curves,  which  are  seen  to  be  good  fits  of  the  numerical 
points .  (Slightly  larger  values  of  t  would  be  expected  for  the  ID/OD  =  1/3 
case  than  for  ID/OD  =  2/3,  but  no  significant  difference  is  apparent  from 
the  Yp  vs  t  plots) . 

The  Vt/D  term  in  eqn.  (2-13)  for  the  hollow  aluminum  rods  impacting 
at  4.6  km/sec  is  therefore  0.95.  By  comparison,  the  Vt/D  term,  for  solid 
rods  at  this  velocity  is  1.5.  Hence  the  hollow  rods  will  be  more  efficient 
by  approximately  0.5  D.  That  is,  a  greater  portion  of  the  length  of  a 
hollow  rod  will  remain  after  the  transient  phase  is  completed  than  will  be 
the  case  for  a  solid  rod  of  the  same  diameter  (or  mass  per  unit  length) .  The 
rod  loss  during  steady  state  penetration  will  be  identical  for  hollow  and 
solid  rods . 


4.3  DISSIMILAR  HOLLOW  ROD  IMPACTS 

One  numerical  solution  considered  the  impact  of  a  hollow  steel  rod 
into  a  T/D  =  2  aluminum  target.  The  results  are  chiefly  of  interest  in  con¬ 
firming  the  applicability  to  dissimilar  material  impacts  of  the  rod  loss 
equations  developed  in  Section  2.2. 

Figure  4-4  shows  the  Yp  vs  t  data  obtained  from  the  numerical  solution. 
The  data  have  been  fitted  with  eqn.  (2-9),  using  values  of  Vco,  the  stagna¬ 
tion  point  velocity,  which  from  eqn.  (2-6)  is 


(4-1) 


Case  23,  ID/OD  =  1/, 


FIGURE  4-4  Penetration  Depth  vs  Time  For  Hollov^ 
Iron  Rod  Impact  into  Aluminum  at 
4.6  km/sec, Target  Thickness/ Rod 
Diameter  =  2. 


impacts  sy 


'  ^T-*T'  — ^  -125X.6-75  (4-2) 
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The  fit  which  is  obtained  using  these  values  in  eqn.  (2-9)  is  seen  to  be 
excellent. 
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APPENDIX  A 

SPATIAL  PLOTS  FROM  NUMERICAL  SOLUTIONS 


The  two-dimensional  hydrodynamic  PIC  code  was  used  for  all  of  the 
impact  solutions  obtained  in  this  investigation.  The  problems  were  run  on 
the  IBM  7094  computer  at  Eglin  APB,  Florida.  For  ease  of  interpretation, 
the  computer  output  can  be  obtained  in  the  form  of  visual  displays  of  spatial 
plots  by  use  of  the  SC4020  plotter.  Three  types  of  spatial  plots  are 
generally  called  for  -  showing  mass  positions,  velocity  fields,  and  pressure 
fields  at  selected  times.  Mass  position  plots  indicate  the  material  dispo¬ 
sition  by  arrangements  of  dots  representing  masses.  Velocity  fields  show 
average  particle  velocities  for  the  mass  points  in  each  cell  by  means  of 
vectors  originating  at  the  cell  centers.  Pressure  fields  are  shown  by 
constant  pressure  contours  (isobars).  In  this  appendix,  the  development  of 
the  impact  solutions  are  illustrated  by  sequences  of  the  mass  position  and 
velocity  field  plots. 

Figures  1-1  and  1-2  of  the  text  summarize  the  impact  cases  which 
have  been  considered  during  the  study.  Spatial  plots  for  cases  1-8  are 
given  in  Appendix  1  of  reference  1.  Up  until  the  time  that  the  initial  shock 
in  a  given  rod  and  velocity  case  reaches  the  rear  surface  of  a  thin  target, 
the  solution  is  identical  to  an  impact  on  any  thicker  target.  Wherever 
possible,  thicker  target  solutions  in  this  investigation  have  hence  been 
"built”  upon  early  stage  solutions  for  thinner  targets.  Thus  in  the  sequence 
of  impact  cases  where  different  thicknesses  of  targets  were  considered 
(l.e,  T/D  =  1/4,  1/2,  1,  2,  3,  5)  the  initial  conditions  for  the  T/D  =  1/2 
solution  were  obtained  from  the  t/d  =  1/4  solution  at  a  point  just  before 
the  shock  in  the  target  had  reached  the  rear  surface  of  the  T/D  =  1/4  plate. 
Similarly  the  T/D  =  1  solution  was  started  using  conditions  from  the 
T/D  =  1/2  problem  just  before  the  shock  reached  the  T/D  =  1/2  rear 
surface. 


Irr.pAct  a!  -i.o  kr"/ier. 
Pod  Did.  =  cO  cm.  Target 
Thickness  ^  180  cm  (T/D  =  3) 


Figure  A1  (facing)  Initial  Conditions 

The  upper  plot  illustrates  the  initial  disposition 
of  the  mass  points  in  the  computational  grid,  while  the 
lower  plot  shows  the  initial  velocity  field  at  the  instant 
of  impact.  In  these  plots  the  left  ordinate  is  the  axis 
of  symmetry  and  the  line,  y  =  0,  corresponds  to  the 
initial  interface  between  rod  and  target  material.  Equi¬ 
valent  solutions  are  obtained  either  by  moving  the  rod 
towards  the  target,  or  by  moving  the  target  towards  the 
rod. 


To  optimize  the' resolution  obtained  in  the  solutions, 
those  portions  of  the  rod  and  target  materials  which  are 
remote  from  the  impact,  and  which  are  hence  not  yet  active 
in  the  impact  process,  are  not  Included  in  the  field  of  view. 
The  rod  itself  is  always  shown  foreshortened,  with  new 
length  being  added  as  it  advances  into  the  target  plate. 

Hydrodynamic  solutions  can  be  linearly  scaled, 
hence  the  only  important  quantity  in  Figure  A1  is  the 
ratio  between  target  thickness  and  rod  radius.  For 
convenience,  these  problems  were  set  up  using  60-cm 
diameter  rods .  The  results  can  be  scaled  to  any  other 
size  by  multiplying  all  lengths  and  times  by  the  same 
factor.  To  convert  to  0.125  in. -diameter  rods,  one 
need  only  multiply  all  lengths  and  times  by  2. 54/(60)  (8)  = 
0.00529166  ....  This  treatment  does  not  alter  the 
thermodynamic  and  mechanicrl  variables  necessary  to 
describe  the  results. 
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IMPACT  CASE  11:  A1  vs  A1  C-  4.  6  km/sec,  T/D  =  j 


Figure  A2  (facing):  Velocity  Fields  at  5.7  and  10.7  usee 


These  plots,  as  well  as  those  in  Figures  A3-A6  ; 

which  follow,  are  taken  from  the  T/D  =  i  ,  i  ,  1,  and  2  * 

solutions  in  Reference  1.  They  are  included  here  to  pro-  ’ 

vide  the  "built-up"  basis  for  the  continuing  solutions  of  _  I 

the  T/D  =  3,5  impacts.  ^ 

The  early  stc^ges  of  the  flow  field  are  seen  in  - 

Figure  A2 .  ^ 

Initially,  most  of  the  field  is  one-dimensional, 
with  shocks  propagating  in  opposite  directions  from  the 
original  interface.  Between  the  diverging  shocks  is  a 
constant  state  region  moving  at  half  the  impact  velocity. 

Pressure  relief,  in  the  form  of  a  rarefaction  fan,  propagates  | 

into  the  high  pressure  region  from  the  free  surfaces.  This  f 

rarefaction  system  eventually  destroys  the  one-dimensional  | 

nature  of  the  initial  flow.  | 
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IMPACT  CASE  U: 


A1  vs  A1  *3  4.0  km/spc ,  T/D  =  3 


Figure  A3  (facing);  Velocity  Fields  at  16.2  and  2  6.2  usee 


These  plots  show  the  continuing  developraent  of  the 
velocity  field.  The  particle  velocities  flowing  radially  into 
the  target  material  interact  with  the  front  surface  to  cause 
upward  flow  — the  front  surface  splash  which  is  experi¬ 
mentally  observed  in  impacts. 
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IMPACT  CASE  11:  A1  vs  A1  @  4.6  km/sec,  T/D 
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IMPACT  CASE  11:  A1  vs  A1  @  4.6  km/sec,  T/D 


IMPACT  CASE  12:  Aluminum  Rod  vs  Aluminum  Plate,  Normal 

Imoact  at  4.6  km/sec.  Rod  Diameter  =  60  cm. 
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IMPACT  CASE  21:  Hollow  Aluminum  Rod  vs 
Aluminum  Plate,  Normal  Impact  at 
4.  6  km/sec,  Rod  OD  -  60  cm, 

ID  t  20  cm  (ID/OD  =  1/3)  ,  Target 
Thickness  =  120  cm  (T/OD  =  2) 


Figure  A1 5  (facing):  Pressure  Fields  at  13.9  and 

18. 6  usee 


The  left  axis  In  these  plots  Is  the  axis  of  symmetry  (the 
center  of  the  hollow  rod).  The  contour  lines  are  10,  50,  100, 
300,  and  500  kb  Isobars.  In  these  plots,  the  pressure  field  is 
expanding  from  the  initial  impact  surface  of  the  rod  wall  on  the 
target.  Convergence  causes  the  pressures  to  be  higher  near 
the  axis  than  in  the  outer  direction.  The  shock  consequently 
propagates  more  rapidly  in  the  inner  direction.  This  is  evident 
at  13.9  psec,  in  that  while  the  pressure  field  extends  nearly  to 
the  axis,  a  distance  of  10  cm  from  th.e  inner  wall  of  the  rod,  it 
has  propagated  only  about  8  cm  outward  from  the  outer  wall  of 
the  rod. 


By  18,6  usee,  the  converging  wave  has  reached  the 
axis,  resulting  in  a  sharp  build-up  of  pressure  there. 
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IMPACT  CASE  22;  Hollow  Aluminum  Rod  vs  Aluminum  Plate, 
Normal  Impact  at  4.6  km/sec. 
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IMPACT  CASE  32:  Aluminum  Plate  Projectile  vs 
Aluminum  Plate  Target, 

30°  Oblique  Impact  at  4.6  km/sec, 
(Planar  Solution) , 

Projectile  Thickness  ("D")  =  60  cm. 
Target  Thickness  =  120  cm 
(T/D  ^  2) 


Figure  A38  (facing):  Initial  Mass  Positions 


Analysis  of  oblique  impact  of  a  rod  requires  use  of 
three  spatial  dimensions.  To  reduce  such  problems  to  two 
dimensions,  a  planar  approximation  (no  z-direction  dependence) 
is  used.  The  solution  therefore  represents  the  edge-on  impact 
of  an  infinitely  long  plate  onto  another  plate.  The  plots  shown 
are  cross-sections.  The  projectile  center  line  is  at  1030  cm. 
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